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Using the plane wave pseudopotential method we performed density functional theory calculations on the stability 
of steps and self-diffusion processes on Ag(lOO). Our calculated step formation energies show that the {111}- 
faceted step is more stable than the {110}-faceted step. In accordance with experimental observations we find 
that the equilibrium island shape should be octagonal very close to a square with predominately {lll}-faceted 
steps. For the (100) surface of fee metals atomic migration proceeds by a hopping or an exchange process. For 
Ag(lOO) we find that adatoms diffuse across fiat surfaces preferentially by hopping. Adatoms approaching the 
close-packed {lll}-faceted step edges descend from the upper terrace to the lower level by an atomic exchange 
with an energy barrier almost identical to the diffusion barrier on flat surface regions. Thus, within our numerical 
accuracy (« ±0.05 eV) there is no additional step-edge barrier to descent. This provides a natural explanation 
for the experimental observations of the smooth two-dimensional growth in homoepitaxy of Ag(lOO). Inspection 
of experimental results of other fee crystal surfaces indicates that our result holds quite generally. 

PACS numbers: 68.55.-a, 68.35.Fx, 68.35. Bs, 71.45. Nt 



I. INTRODUCTION 

The study of the morphology of surfaces is one of the 
oldest areas of crystal growth and has attracted particu- 
lar attention in recent years due to advances in epitaxy 
and the possibility of controlled growth of new materials 
and bulk quality surfaces. The structure of an epitax- 
ial film is often thought to be determined by thermody- 
namic principles, in which the film is assumed to adopit 
the energetically stable thermal-equilibrium structure.EJ 
However, there are many examples which demonstrate 
that under typical growth conditions growth is often gov- 
erned by kinetics rather than thermodynamics. A lot of 
experimental works show that the growth mode (the evo- 
lution of film structure with coverage) can be altered by 
varying the growth conditions such as the substrate tem- 
perature and deposition rate, or by introducing defects or 
surfactants.H A notable example is the study of homoepi- 
taxial growth of Ag(lll): Under clean surface conditions 
the growing surfaces are very rough with mormtains as 
high as 30-40 atomic layers (multilayer growth) .I3'D These 
studies revealed that the film structure is governed by the 
kinetics of surface diffusion. Hence, knowledge of the ki- 
netics such as mass transport across terraces and steps 
and the diffusion parallel to steps is essential to reach 
understanding of the morphology and quality of growth. 

To date most calculations of atomic diffusion on sur- 
faces employ computationally fast semiempirical meth- 
ods. They have been widely used by methodological 
simplicity and under the assumption that the results, 
although quantitatively not very accurate, may still ex- 
plain general trends. The main approximation of these 
methods concerns the kinetic energy operator of the elec- 
tron and the neglect of a self consistent adjustment of 
the electron density. Thus the essence of the mechanism 
governing the breaking and making of chemical bonds is 



missing. Obviously, quite often, as e.g. in the example 
discussed in this paper, quantitative results determine 
the qualitative features. For example, when the energy- 
barrier at the step edge is bigger than the diffusion bar- 
rier at the flat surface the growth mode is rough, and 
when the step-edge barrier is equal or smaller, the growth 
mode is flat. The proper treatment which gives quantita- 
tively and qualitatively more accurate results and takes 
the quantum mechanical properties and the chemistry of 
interatomic interactions properly into account is density 
functional theory (DFT). 

In this paper we report DFT total energy calculations 
of step formation and diffusion of Ag adatoms on the 
(100) surface of fee silver, which extends work presented 
earlier.u 

Homoepitaxy of silver exhibits a different behaviour 
in the growth mode depending on substrate orienta- 
tions: The growing surface of Ag(lOO) is smooth and flat 
whereas that of Ag(lll) is typically rather rough. For 
Ag(lOO) surfaces reflection high-energy electron diffrac^. 
tion (RHEED) intensity oscillations have been observed^ 
in the temperature range from 200 to 480 K, which 
indicates that the growth proceeds by a smooth two- 
dimensional (2-D) mode. In contrast, for Ag(lll) sur- 
faces RHEED intensity oscillations are absent at all 
temperatures, suggesting that there is no 2-D growth 
but that the surface becomes very rough upon growth. 
This finding for Ag(lll) was confirmed by x-ray refiec- 
tion experimentaJ and. -scanning tunneling microscopy 
(STM) measurementsJjQ The former study shows that 
the growth mode changes from multilayer to step-fiow 
growth over the temperature range from 175 to 575 K. 
The latter studies show an STM topography of "moun- 
tain and canyon" -like structures with heights of 30-40 
atomic layers. 

^From thei*. analysis of STM images for Ag(lll) Vri- 
jmoeth et aln concluded that the multilayer growth is 
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attributed to an additional step-edge barrier. This is 
the barrier for a diffusing adatom descending steps mi- 
nus the surface diffusion barrier. The authors obtained 
that ihe additional energy barrier is 0.15 eV. Bromann 
et alB developed a method to determine with high accu- 
racy attempt frequencies and activation energies for ter- 
race diffusion as well as for the step-down process. This 
analysis was done by measuring the nucleation rate on 
top of islands as a function of island size and tempera- 
ture. They obtained for the additional energy barrier a 
value of 0.12 eV very close to the 0.15 eV by Vrijmoeth 
et al. Recent STM studies by Morgenstern et alB find 
the additional step-edge barrier to be 0.10 eV. This ad- 
ditional energy barrier which adatoms encounter at step 
edges explains why the growth of silver perpendicular to 
the (111) surface proceeds by a very rough multilayer 
mode. Interestingly, as noted above, the growth of the 
(100) surface is qualitatively different. As the calcula- 
tions in Rcf. ^ showed the additional step-edge barrier 
is practically zero. Below we extend our previous study 
and repeat convergence tests. 

Furthermore, we consider the shape of islands formed 
on Ag (100) in thermal equilibrium. There are two types 
of monolayer- high steps for fee (100) surfaces as shown 
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FIG. 1. Top view of monolayer-high steps on a fee (100) 
terrace: a {lll}-microfaceted step running along the (110) 
direction and a {110}-microfaceted step running along the 
{ 100 ) direction. The two step edges are labeled as 
(110) /{111} and (100) /{HO}, respectively. For our ex- 
ample Ag(lOO) the nearest-neighbor spacing, i.e., a surface 
lattice constant as, is 2.92 A. 



m Fig. |1|: One is a close-packed step running along the 
(110) direction which has a {111} microfacet at edges 
and the other an open step running along the ( 100 ) di- 
rection which has a {110} microfacet at edges. A ratio 
of formation energies of these two steps determines the 
equilibrium shape of islands on Ag(lOO) (as obtained by 
the Wulff construction) B Here we apply the ab initio to- 
tal energy calculations to the energies of these two sorts 
of steps. 

The paper is organized as follows. In Sec. II, we de- 
scribe the method of calculation. In Sec. Ill, the equilib- 
rium shape of islands is discussed as it follows from our 
calculated step-formation energies. In Sec. IV we present 
results for diffusion processes of a Ag adatom on flat and 
stepped Ag(lOO) surfaces, and discuss their implication 
on epitaxial growth. Finally the paper is concluded by 
Sec. V. 



II. METHOD OF CALCULATION 

The total energies and forces on atoms are calculated 
by the DFT. We use a norm-conserving pseudopotential 
method together with the Car-Parrinello-like technique 
for iterative minimization of the Kohn-Sham total en- 
ergy.Ej The electronic wavefunctions are expanded in a 
plane wave basis set. Most calculations have been car- 
ried out primarily using the local-density approximation 
(LDA) for the exchange-correlation functional.llil At im- 
portant geometries we repeat the calculations using the 
ge:^ralized gradient approximation (GGA) of Perdew et 

We use a repeating slab geometry to simulate the ac- 
tual surface. A (3x3) periodicity is employed in the 
lateral directions, which we tested gives that artificial 
adatom-adatom interaction is sufficiently weak and in 
fact negligible for the questions of concern. An artifi- 
cial periodicity is imposed along the perpendicular di- 
rection to the surface to construct a three-dimensional 
unit cell. We take a slab of three atomic layers and of 
a 10.35 A vacuum region. This rather thin slab is afe 
ceptable, as we adsorb the adatom only on one side.Ej 
We relax the adatom and the top-layer atoms according 
to a damped Newton dynamics, while the other atoms 
in the two bottom layers are kept at the bulk positions. 
The optimized geometries are identified by the require- 
ment that the remaining forces acting on the atoms are 
smaller than 0.05 eV/A. 

The LDA (GGA) ionic pseudopotentials including the 
tightly bound 4d-electron state as valence states are cre- 
ated using the nonrelativis.tic (semirelativistic) scheme 
of TrouUier and Martins. Eil The fully separable norm- 
conseapj|ig_ pseudopotentials of the Kleinman-Bylander 
fornxJO'ED are constructed with the s pseudopotential 
as the local component. We use a kinetic energy cutoff 
of 40 Ry to expand the wavefunctions in plane waves. 
The k-space integration is performed using nine equidis- 
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TABLE I. Structural properties of fee silver obtained 
from the plane wave pseudopotential and all-electron lin- 
earized-augmented-plane-wave (LAPW) methods, for nonrel- 
ativistic (NR) and semirelativistic (SR) calculations in the 
LDA and GGA. The T —+ K experimental data are given 
for comparison. 
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(SR) 


(SR) 


(SR) 




ao (A) 


4.14 


4.18 


4.00 


4.17 


4.07 


Bo (GPa) 


99 


90 


136 


85 


102 



^Reference |l| 

'^The experimental values are taken from Reference lli 



tant k points in the surface Brillouin zone (SBZ) of the 
(3x3) cell. To improve the k-space integration and sta- 
bilize convergence the electronic states are occupied af-, 
cording to a Fermi distribution with k^Td =0.1 eV,t2l 
and then the total energies are extrapolated to zero elec- 
tronic temperature. In order to attain the iterative so- 
lutions of the Schrodinger equation we start with initial 
wavefunctions obtained from the self-consistent solutions 
to the Kohn-Sham Hamiltonian in a mixed basis set of 
pseudo atomic orbitals and plane waveSrwith a cutoff en- 
ergy of 4 Ry, as developed by Kley et alS3 This approach 
gives much faster convergence than the conventional one 
starting with the initial wavefunctions created from the 
diagonalization of the Hamiltonian and a small cutoff en- 
ergy in the plane wave basis set. 

To test the accuracy of the pseudopotentials, we cal- 
culated the ground state properties of bulk silver, using 
512 k points in the Brillouin zone. The results derived 
from the Murnaghan equation-of-state fit to the calcu- 
lated data are listed in Table |. Zero-point vibrations 
are not considered in these theoretical results. The cal- 
culated values are in good agreement with T — > K 
experimental data. 

III. STEP FORMATION ON (100) SURFACES 

A. Structure and energetics of the flat surface 

Before we describe the step formation on Ag(lOO) it 
is worthwhile to study the properties of the flat Ag(lOO) 
surface. The surface formation energy per surface atom 
is given as 

= (-Eslab - NAg^J,Ag)/Nsu,:^ , (1) 

where i?siab is the total energy of the slab of A^Ag silver 
atoms and iVgurf is the number of atoms at the surface. In 
thermal equilibrium the silver chemical potential is equal 
to the energy of a silver atom in the bulk. In order to 
achieve an optimum of error cancellations the calculation 
of the bulk energy /iAg is done with the same k-points 
sampling as in the slab calculations. The work function 



TABLE II. Calculated results for the properties of flat 
Ag(lOO): surface formation energy a, work function (f), and 
multilayer relaxation Adij relative to the interlayer spacing in 
the bulk do- It is noted that the experimental surface energy 
is an estimate from liquid surface tension measurements. 
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(%) 


This work (NR-LDA) 


0.57 


4.39 


-2.2 


0.4 


FP-LMTO (NR-LDA)'' 
FP-LMTO (SR-LDA)'' 


0.63 


4.43 


-1.9 




0.62 


4.77 


-1.9 
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0.65" 


4.64'^ 


0± 1.5° 


0± 1.5'= 



^Reference 20| 
''Reference p 
''Reference W 
''Reference 23^ 
"Reference 24\ 



is given by the difference between the electrostatic energy 
in the middle of the vacuum region and the Fermi energy. 
The multilayer relaxation Adij is defined by the change in 
spacing between layers i and j compared to the interlayer 
spacing in the bulk, do- 

We calculate the properties of flat Ag(lOO), using the 
slab of 5 atomic layers, a vacuum region equivalent to 5 
atomic layers, and 96 k points in the SBZ of the (1x1) 
cell. In this calculation we used the lattice constant 
4.14 A of k-sampling-converged calculation which dif- 
fers only by 0.1 % from the calculated lattice constant 
using the above k-points. The results are sumraarized 
in Table The results in our previous papeitl were 
obtained using a slab of 4 atomic layers, a vacuum re- 
gion equivalent to 4 layers, and 64 k-points in the SBZ 
of the (1x1) cell. Due to these changes the present 
results differ slightly from those of Ref. ||. The calcu- 
lated results are in good agreement with a previous full- 
potepti^d linear-muffin-tin-orbital (FP-LMTO) calcula- 
tion.Elt2l Our calculations of the lattice relaxation give 
Adi2/do = -2.2% and Ad2z/da = 0.f%. The low-energy 
electron diffraction (LEED) analysis p3 in which the first 
(di2) and second (^23) interlayer spacings were deter- 
mined, shows small relaxations with A(ii2/do = 0± 1.5% 
and A(i23/do = 0±1.5%. This result is in good agreement 
with the present DFT-LDA results within the experimen- 
tal accuracy, (see also the. influence of surface vibrations, 
which is neglected here£3). 

B. Structure and energetics of steps 

In this subsection we consider the stability of the two 
most densely packed monolayer-high steps on Ag(lOO): 
the close-packed (110)/{111} and open (100)/{110} 
steps (see Fig. |^). To compute step energies we use the 
slab geometry of the grooved surface sketched in Fig. ||. 

The simulation cell has the periodicity of along the 
step edge and of Ly along the perpendicular direction to 
the step edge. We take to be Os and \/2as for {111}- 
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FIG. 2. Schematic representation of the A^-layer slab of the 
grooved surface. Periodic boundary conditions are applied 
along the a:: and y axes. 



and {110}-faceted steps, respectively. Ly is treated as a 
variable, so that the ledge separation / is equal to Ly/2 
(see Fig. ^). The formation energy of the grooved surface 
E'surf is extracted from the calculations of the A^-layer 
slab shown in Fig. |[ by applying the formula 

£;surf(A^) =i^slab(iV)-iVAg/^Ag • (2) 

The step formation energy per unit length is 

\{N)^{E,^,i{N)-N,^,ia{N))lL, . (3) 

In Eq. ^, (t{N) is the surface energy per surface atom 
from the calculations of the A^-layer slab of the flat (100) 
surface and iVgurf is the number of surface atoms at the 
grooved surface: Nguri = L^Ly/a^. The surface energy 
is evaluated by averaging the surface energies of the thin 
and thick region of the slab, i.e. aa,vc{N) = (1 — a)cr(A^) + 
a(T{N + 1) (a is a ratio of the number of surface atoms in 
the thick region and the total number of surface atoms). 
The averaging improves the accuracy of the calculated 
step energies. 

Several calculations were done to test the dependence 
of the calculated step energies on the ledge separation 
I and the number of atomic layers N. The results of 
the {lll}-faceted step are summarized in Table [IL It 
shows that the step energy decreases with increasing the 
ledge separation I. In particular, the step energji—X can 
be expanded with respect to the ledge separations 



A = Ao + A/f 



(4) 



In this equation the first term is the energy of an isolated 
step, and the last term gives the effective interaction be- 
tween steps. This term includes possible energetic contri- 
butions such as dipole-dipole or elastic interactions. For 
the 3-layer grooved surface we obtain Aq — 0.150 eV/a^ 
and A = 0.041 eV/a^ in the relaxed geometry by fitting 
the variation of the step energy with the ledge separation 
I to Eq. H We find that the step energy already reach the 
isolated one at I = Sos with an error of 0.005 eV/ag. The 
step energy decreases by 0.024 eV/as with increasing the 



TABLE III. Convergence of the formation energy of the 
{lll}-faceted step in the relaxed geometry as a function of 
ledge separation I and the number of atomic layers A'^, calcu- 
lated by DFT-LDA. 
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E.urf (eV) [eq§ 


4.951 


7.258 


9.593 


7.172 


o"avc (eV/atoni) 


0.579 


0.579 


0.581 


0.576 


A (eV/a.) [eq| 


0.160 


0.155 


0.152 


0.130 



slab thickness from 3 to 4 layers. We find that, while 
the step energy changed by 20 % as the slab thickness 
is increased from 3 to 4 layers, a ratio of step energies 
changes only by 5 %. Thus, we see that the slab of 4 
layers is sufficient to obtain reasonably converged result 
of a ratio of step energies. 

Using the 4-layer slab of the grooved surface and the 
ledge separation of 3as, we obtain the formation energy 
of 0.130 eV/fls (LDA) for the close-packed {lll}-faceted 
step. For the {110}-faceted step we take a 4-layer slab 
and ledge separation of 2\f2as. The obtained step en- 
ergy is 0.156 eV /as (LDA). Thus, the close-packed {111}- 
faceted step is more stable than the open {110}-faceted 
one. This confirms the experimental observation that 
the {lll}-faceted steps are preferentially formed in ther- 
mal equilibrium. This result is not unexpected because 
of atomic geometries, i.e. the local coordinations of 
the two steps differ noticeably: Step-edge atom at the 
{lll}-faceted step has seven nearest neighbors, while 
the step-edge atom at the {110}-faceted step has ouli' 
six neighbors. Using the AFW ( Adams- Foiles-Wolfer^ll 
EAM (embedded atom method) functions. Nelson et aZ.E2l 
found A to be 0.102 eV/a^ and 0.142 eV/os for {111}- and 
{llOWaceted steps, respectively. With the VC (Voter- 
ChenO) functions they obtained A to be 0.135 eV/os for 
{lll}-faceted step. 

When the formation energies of all steps were known, 
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FIG. 3. Calculated shape of islands in thermal equilibrium. 
The close-packed {lllj-faceted edges are longer with an edge 
length ratio l'^^oo)/{iiq} . ^(iio>/{iii} 3.;^q ^j^^^ ^^le open 

{110}-faceted ones. 



4 



the equihbrium shape of an island can be obtained by 
finding the island shape with miaimum free energy, or by 
applying the Wulff constructions For the (100) surface it 
is certainly plausible that only two steps are important. 
Thus we expect octagonally shaped islands as shown in 
Fig. H The calculated step formation energies imply that 
the {lll}-faceted edges should dominate and that the 
edge-length ratio should be l(ioo>/{iio} . 2,(iio)/{iii} ^ 

3 : 10. This theoretical finding can be compared to exper- 
imental observations. In fact, for ion-bombarded Ag(lOO) 
surface the preferential formation of the close-packed 
{lllj-faceted step edges has been observed at room tem- 
perature,E2l and for Ir(lOO) Chen and Tsong, using the 
field ion microscope, also found that the equilibriuia is- 
land shape is rectangular with {lll}-faceted steps.E2l 

Next we consider the relaxations of the stepped sur- 
face. Figure Q shows the geometry optimized structure 
of {lll}-faceted steps. The relaxations of surface atoms 
in the upper terrace are nearly identical to those of sur- 
face atoms of flat Ag(lOO): surface atoms 1, 1', and 1" 
relax toward the bulk by ^ 2% of the bulk interlayer 
spacing c?o [see Fig. ^ (a)]. Atoms 2 and 2' at step bot- 
tom slightly relax upward and toward atoms 3 and 3', 
but these relaxations of lower terrace atoms are nearly 
negligible. Figure || (b) displays the bond lengths of step- 
edge atoms with their underlying nearest neighbors and 
the numbers show that the bond lengths are very close to 
those between surface atoms and their nearest neighbors 
in the subsurface sites of flat Ag(lOO). Similarly, at the 
grooved surface with the {110}-faceted steps downward 
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FIG. 4. (a) Side and (b) top views of stable geometry of the 
grooved surface of the ledge separation Sos with {lll}-faceted 
steps. In (a) arrows and numbers in unit of A indicate di- 
rections and magnitudes of atomic displacements from the 
unrelaxed ideal step geometry. The percentage changes in 
distances to nearest-neighbors of step-edge atoms relative to 
the nearest-neighbor separation in bulk are shown in (b). 



relaxations of upper terrace atoms dominate. 

IV. SELF-DIFFUSION ON FLAT AND STEPPED 
(100) SURFACES 

A. Adsorption and self-diffusion on flat surfaces 

We now present the results of adsorption and diffusion 
of a Ag adatom on the flat Ag (100) surface. For met- 
als, metal adatoms on surfaces preferentially adsorb at 
the highest coordination site. Indeed the equilibrium ad- 
sorption site for a Ag adatom on Ag (100) is the fourfold 
hollow [Fig. ||(a)]. In the optimized geometry, the four 
nearest neighbors of the adatom slightly distort in the 
lateral directions, opening the hollow site even further. 
The lateral displacement from the ideal clean surface site 
is 0.02 A. The adatom is located 1.84 A above the sur- 
face layer, and the bond length between the adatom and 
its neighbors is 2.78 A, i.e., 5 % shorter than the in- 
teratomic distance in the bulk. This follows the typi- 
cal trend, namely that bond strength per bond decreases 
with coordination and correspondingly, bond length in- 
creases with coordination. 

In metal-on-metal surface diffusion, an adatom typi- 
cally moves across a flat surface by a series of hops be- 
tween adjacent equilibrium adsorption sites. For fee (100) 
an adatom in a fourfold site moves over the twofold 
bridge site into a neighboring fourfold equilibrium site 
[see Fig. ^(b)]. An adatom may visit all of the fourfold 
surface sites, forming a (1 x 1) square pattern. An al- 
ternative mechanism for surface diffusion is atomic ex- 
change where a diffusing atom moves by displacing a 
neighboring surface atom. For the fee (100) surface the 
exchange process of an adatom with a surface atom oc- 
curs preferentially in the [010] and [001] directions. At 
the saddle point of the transition the adatom and sur- 
face atom are above a vacated surface-layer site [see 
Fig. 11(c)] . Subsequently the displaced surface atom be- 
comes a new adatom in a next-nearest-neighbor fourfold 
site. When this diffusion mechanism is active, an adatom 
may visit only half of the fourfold surface sites, forming 




[Oil] 

FIG. 5. Top view of an adatom at (a) the fourfold hollow 
site, (b) the transition state for the hopping diffusion (the 
twofold bridge site), and (c) the transition state for the ex- 
change diffusion. The white and gray circles represent the ad- 
and surface atoms, respectively. 
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a c(2 X 2) square pattern. Experimental evidences of the 
self-diffusion by exchange displacement on fee (100) haj^e 
been obtained by Kellogg and Feibelpan on Pt (100)e3 
and by Chen and Tsong on Ir (100)Ej. The process haa 
been theoretically studied by Feibelman for Al (100). til 
His total-energy calculations of the self-diffusion barrier 
on Al (100) show that atomic diffusion via exchange is en- 
ergetically favored by a factor of three over hopping. It is 
an open question if the exchange diffusion plays a role for 
other fee tiJQO) surfaces as well. The EAM calculations of 
Liu et alM suggest that, besides Pt(lOO) and Ir(lOO), 
exchange displacements are also favored for self-diffusion 
on the (100) surfaces of Pd and Au. Fcp-^elf-diffusion on 
Cu(lOO) effective medium calculationg2a predicted that 
exchange displacement occurs preferentially, but this re- 
sult is in contrast to other_aemiempirical studie£3 and 
to a DFT-LDA calculationO 

Our total-energy calculations for the self-diffusion on 
Ag (100) clearly show that the hopping mechanism is en- 
ergetically favored. The energy barriers for hopping diffu- 
sion are 0.52 eV (LDA) and 0.45 eV (GGA) which, is close 
to the result of the FP-LMTO LDA calculationEd , which 
got 0|^ eV. Using the AFW and VC EAM functions Liu 
et alSB obtained a similar value, namely 0.48 eV. Our 
results show that the use of the GGA functional for the 
exchange-correlation energy lowers the energy barrier by 
13 % (0.07 eV) compared to the LDA resuh. Recent 
STM studies of homoepitaxy on Ag(100)Ea, combined 
with Monte-Carlo simulations of an appropriate model 
for fee (100) homoepitaxy, determined the diffusion bar- 
rier to be 0.33 eV. Low-eaergy ion scattering measure- 
ments of Langelaar et alx^ obtain a value of 0.40 eV. 
At this point we note that the experimental analyses of 
the diffusion barriers are not without problem. Thus, 
we conclude that our GGA result {Ed = 0.45 eV)B is 
well in agreement with (subsequently obtained) available 
experimental information. In the optimized geometry of 
the bridge site, the saddle point for the hopping diffusion, 
the two nearest neighbors of the adatom are pushed away 
and downwards by 0.05 A. The adatom is located 2.33 A 
above the surface layer. The bond length between the 
adatom and its two neighbors is 2.69 A, i.e. 3 % shorter 
than in the fourfold hollow. This again follows the well 
known trend: Each of the two bonds at the bridge site is 
stronger than each of the four bonds at the hollow site. 

To obtain the energy barrier for the exchange diffusion 
we compute the total energy of the transition state shown 
in Fig. ||(c). We obtain an energy barrier of 0.93 eV 
(LDA) and 0.73 eV (GGA) which is much higher than 
that for the hopping process. Thus, we see that Ag 
adatoms diffuse across flat regions of Ag (100) by hop- 
ping. It is interesting that for the exchange geometry the 
difference between LDA and GGA is noticeable (0.20 eV). 
The GGA value is comparable to the EAM results. Us- 
ing the AFW and VC EAM functions Liu et al. found 
an energy barrier of 0.75 eV and 0.60 eV for the ex- 
change process, respectively. In the optimized geome- 
try of the transition state for the exchange diffusion [see 



rfj2= 2.75 A 



(a) 



(b) 



(c) 




FIG. 6. Valence charge density at the transition state for 
the exchange diffusion. A schematic plot of the optimized 
atomic geometry is shown in (a). Panels (b) and (c) display 
the electron density. In (b) the plot is in the plane cutting 
atoms 5, 1, 2, and 6. In (c) the plot is in the plane cutting 
atoms 3, 1, 2, and 4. The value of the lowest-density contour is 
3xl0~^ e/(bohr)^ and subsequent contours differ by a factor 
of 2.2. Solid circles mark the second layer atoms. 

Fig. ^(a)], each of the two top- most adatoms has four 
bonds with neighbors. The distance between the two top- 
most adatoms is 2.75 A, 1 % shorter than in the fourfold 
hollow. These two atoms are located 1.52 A above the 
surface layer. The four atoms neighboring these atoms 
are pushed away by 0.06 A and downwards by 0.05 A. 
The bond length between the top-most and underlying 
Ag atoms is 2.71 A, 3 % shorter than in the fourfold hol- 
low. Another bond length between the atoms 1 and 5 
(between the atoms 2 and 6) in Fig. |(a) is 3.07 A, 10 % 
longer than in the fourfold hollow, and 5 % longer than 
in the bulk. We note that at the (100) surface of alu- 
minum, the exchange displacement is supported by the 
formation of covalent bonds, which is made possible (and 
indeed plausible) by the sp valence of the group III el- 
ement Al. For silver we find that the electron density 
does not reflect a pronounced covalency effect as shown 
in Fig. ^(b) and (c), which explains why the exchange 
displacement does not play a role at flat Ag(lOO). 



B. Adsorption and self-diffusion on stepped surfaces 

This subsection is concerned with the diffusion of a 
Ag adatom at the {lll}-faceted steps which are prefer- 
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number of atomic rows within the (lOO)-oriented terraces is 

three. Distance(aJ 



entially formed on Ag (100) surfaces. The calculations 
were done for a high index (511) surface which is vicinal 
to (100) (see Fig. 0). This surface consists of (100) ter- 
races which are three atom rows wide. The step edges 
are perpendicular to the [100] and the [Oil] directions. 
The periodicity along the step edge is taken to be three 
surface lattice constants. Using a vicinal surface allows 
to consider a smaller number of atoms within the unit 
cell than that required for the ordered step arrays on 
flat surfaces (the grooved surfaces). Using the vicinal 
(511) surface of 4 atomic layers, we obtain the formation 
energy of 0.136 eV/os (LDA) for the close-packed {111}- 
faceted step, which is close to the corresponding value 
0.130 eV/os (LDA) from the 4-layer slab calculations of 
grooved surface (see Table III). In Fig. ||(b) we display 
the total energy of a Ag adatom diffusing on the stepped 
Ag (100) by the hopping process. For any position of the 
path indicated in Fig. ||(a) the total energy is obtained 
by relaxing all coordinates of the atoms of the surface 
layer as well as the height of the diffusing atom. 

Figure ||(b) clearly shows two stable adsorption sites: 
The hollow site (M) at the step, at which the adatom 
is fivefold coordinated and the hollow site (H) on the 
terrace, at which the adatom coordination is four. Due 
to the higher coordination it is indeed plausible that the 
DFT results give that the adsorption site AI is more sta- 
ble than site H. The energy difference between the M 
and H sites is 0.43 eV (LDA) and 0.32 eV (GGA). In the 
optimized geometry of the site M the distance between 
the step-edge atom and adatom is 2.89 A, 1 % shorter 
than in the bulk. The bond length between the two step- 
bottom atoms and adatom is 2.78 A, and that between 



FIG. 8. Total energy of a Ag adatom diffusing along the 
indicated path by a hopping (roll-over) process, calculated 
within the LDA. A top view of the vicinal (511) surface is 
shown in (a). The {lll}-faceted edges are aligned along the 
[Oil] direction. The distance is given in unit of the surface 
lattice constant as. 



the two lower-terrace atoms and adatom is 2.84 A. It is 
also interesting to note that desorption of an adatom at 
the step edge to the flat region requires to overcome a 
very high energy barrier (LDA: 0.96 eV; GGA: 0.76 eV). 

When the adatom rolls over the ledge from H site on 
the upper terrace to M site on the lower terrace, we ob- 
tain an energy barrier of 0.70 eV (LDA) and 0.55 eV 
(GGA). Thus, we find the "additional step-edge barrier" 
to be 0.18 eV (LDA) and 0.10 cV (GGA) to step down 
from the upper terrace by the hopping (roll-over) process. 
This value is by about 30 % (20 %) higher in the LDA 
(GGA) than the diffusion barrier at the flat region. At 
the transition state for the roll-over process the adatom is 
located near the bridge site at the ledge [Sh in Fig. ||(a)]. 
In the optimized geometry the two nearest neighbors of 
the adatom are slightly pushed away, and the adatom is 
located 2.27 A above the upper terrace. The bond length 
between the step-edge atoms and adatom is found to be 
2.72 A, 1 % longer than in the bridge site at the fiat re- 
gion and 7 % shorter than that of a bulk atom: Each 
of the two bonds at site Sh is weaker than that at the 
bridge site on the fiat region. At the bridge site on the 
fiat terraces the adatom has four next-nearest-neighbors, 
while at Sh the adatom has two next-nearest- neighbors. 
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FIG. 9. Total energy of a Ag adatom diffusing across a step 
by an exchange process, calculated within the LDA. In (a) the 
adatom 1 is adsorbed at the fourfold hollow site H. The total 
energy as a function of the distance of the step-edge atom 2 
from the undistorted step edge is shown in (b). 

Thus, we could understand the additional energy barrier 
in terms of weaker bonds between the adatom and under- 
lying surface atoms, and decrease of the number of the 
next-nearest-neighbors compared to the flat surfaces. 

The other possibility to step down from the upper ter- 
race is via an exchange, where adatom 1 displaces step- 
edge atom 2 and subsequently the latter moves to the 
fivefold coordinated site M [see Fig. 1(a)] . Figure |(b) 
displays the results of the step-down diffusion by the ex- 
change process. Within our numerical accuracy for en- 
ergy differences ±0.05 eV) the calculated energy bar- 
rier of 0.52 eV (LDA) and 0.45 eV (GGA) is almost iden- 
tical to that of the hopping diffusion at the flat region. 
The transition state for the exchange process is identifled 
to be near the bridge site formed by two step-bottom 
atoms 3 and 4 on the lower terrace [Sex in Fig. ^(a)]. 
The geometry of the optimized structure of the transition 
state is given in Fig. ^ The top-most adatom (no. 1) 
is located 1.26 A above the upper-terrace atoms. The 
height is much lower than the value of 1.52 A at the flat 
region. Inspection of the geometry at Sex shows that the 
local coordination of atoms 1 and 2 concerned with the 
exchange displacement remains high, since the atom 1 
follows in close contact to atom 2 during the process: At 
Sex each of the two atoms 1 and 2 has five bonds with 




t/l2=2.74 A 
iil8= 2.69 A d2A= 2.72 A 
d\i= 2.76 A d23 = 2.79 A 
d/l7=2.96A d26= 2.99 A 
dlS=3.40A d25 = 3.45 A 



FIG. 10. Schematic view of the optimized atomic structure 
at the transition state Sex of the step-down diffusion by the 
exchange process, calculated within the LDA. dis denotes the 
bond length between the adatom 1 and second-layer atom 
below. 

neighbors, while at the saddle point of the exchange dif- 
fusion on the flat region each of the two top-most atoms 
has only four bonds [Fig. ||(c)]. The origin of the lower 
diffusion barrier of the step-down motion by the exchange 
displacement compared to the flat region is thus the ad- 
ditional bonds formed at the obtained saddle point. Our 
finding that there is no additional energy barrier to de- 
scend from the upper to the lower terrace provides a nat- 
ural explanation for the smooth 2-D growth of silver in 
the (100) orientation. By analyzing thCpS-TM images of 
Ag(lOO) during deposition Zhang et alEB gave an esti- 
mate of 0.025 ± 0.005 eV for the additional step-edge 
barrier. Within our numerical accuracy (« ±0.05 eV) 
this is in good agreement with our result that the addi- 
tional step-edge barrier is negligible. 

Our finding that for Ag(lOO) an adatom goes down at 
step edges by an exchange process is expected to apply 
for other noble metals and other fee transition metals, 
as well. In fact, for hompepitaxial growth of Cu(100),e2I 
Pd(100),Eland Fe(100)El damped oscillations of the scat- 
tering intensity have been observed. In order to ex- 
plain the intensity oscillations associated with smooth 
2-D layer-by-layer growth, the step-edge barrier for the 
descent of adatoms was assumed to be negligible. For 
a hopping (roll-over) process we expect that an adatom 
generally encounters an additional (positive) step-edge 
barrier [~ 0.1 eV for Ag(lOO)] because of the low coordi- 
nation of the transition state. Thus, we expect that quite 
generally for fee (100) metal surfaces step-down motion 
proceeds by an exchange process, similarly to Ag(lOO). 

Finally we address the atomic motion along the close- 
packed {lll}-faceted step to get ideas about the possible 
roughness of steps and kinetic growth shapes of islands. 
We find that the diffusion along step edge proceeds by 
a hopping process. An adatom in a flvefold coordinated 
adsorption site moves over a bridge site E into a neigh- 
boring fivefold coordinated site [see Fig ^(a)]. ^From 
Fig. |(b) we see that the barrier (LDA: 0.30 eV; GGA: 
0.27 eV) is significantly lower than the surface diffusion 
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barrier Ed (LDA: 0.52 eV; GGA: 0.45 eV). This lower en- 
ergy barrier indicates that Ag islands formed on Ag (100) 
should be compact. Atoms which reach the step edges 
will certainly be able to diffuse parallel to the steps and 
thus local thermal equilibrium is attained. We therefore 
expect rather straight step edges and no fractally shaped 
islands. In the optimized geometry of site E the distance 
between the two step-edge atoms and adatom is 2.99 A, 
2 % longer than in the bulk. We find the bond lengths of 
2.66 A between the step-bottom atom and adatom and of 
2.84 A between the lower terrace atom and adatom and 
note that the coordination varies from five to four as the 
adatom moves from site M' to site E [Fig. ^(a)], while 
for diffusion on flat terraces the coordination varies from 
four to two. Thus, the lower barrier for diffusion parallel 
to steps can be understood in terms of smaller variation 
in coordination. 



V. CONCLUSION 

In the above Section, we presented results of first- 
principles total-energy calculations for the electronic 
structure and energies of steps and for various micro- 
scopic self-diffusion processes at Ag(lOO). 

For fee (100) surfaces two types of monolayer-high 
steps are of particular importance: the close-packed 
{111}- and more open {110}-faccted steps. Our total- 
energy calculations show that the {lll}-faceted step is 
more stable than the {110}-faceted one. This is un- 
derstood in terms of different coordination of step-edge 
atoms. In accordance with experimental observations, 
we find that in thermal equilibrium the shape of island 
should be octagonal but very close to a square with 
{lll}-faceted edges. The edge-length ratio of the oc- 
tagon is calculated as 10:3. The energy barrier for self- 
diffusion of a Ag adatom along {lll}-faceted edges is 
found to be significantly lower than the surface diffusion 
barrier. Thus, we expect rather straight step edges and 
no fractally shaped islands. 

At fiat regions of Ag (100) Ag adatoms are found to dif- 
fuse by a hopping process. The obtained energy barrier 
is 0.45 eV in the GGA which is veryj-close to the ex- 
perimental estimates which are 0.40 eVEJ and 0.33 eVL3. 
In contrast to atomic motion on fiat region, the descent 
of adatoms at steps proceeds by an exchange process. 
The calculated energy barrier is almost identical to the 
barrier at flat region. This indicates that there is no 
additional step-edge barrier to diffuse across step edges. 
This finding is in sharp contrast ta the additional en- 
ergy barrier (A£;,'^f^"^^ = 0.15 evi and 0.12 evi and 
0.10 eVij) experimentally found at Ag(lll) actuating a 
rough growth morphology. The calculated result implies 
good interlayer mass transport at Ag (100) and thus ex- 
plains the smooth 2-D growth experimentally observed 
in homoepitaxy of Ag (100). Inspection of experimental 
works indicates that the step-down motion at steps by 



an atomic exchange takes place quite generally for other 
fee metals. 
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